Defect-induced supersolidity with soft-core Bosons 
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We determine the zero-temperature phase diagram of two-dimensional Bosons with finite-range 
soft-core interactions. For high particle densities, we demonstrate that the ground state can be a 
commensurate density- wave- type supersolid. For low densities, the system is shown to form a solid 
in which superfluidity is provided by delocalized zero-point defects. This provides the first example 
of continuous-space supersolidity consistent with the Andreev-Lifshitz-Chester scenario. 
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Spontaneous symmetry breaking is a focal principle of 
condensed matter physics - yet, simultaneous breaking of 
fundamentally different symmetries represents a rare and 
exotic phenomenon. A prime example is the so-called su- 
persolid state [1 , which displays both crystalline and su- 
perfiuid properties, i.e. the simultaneous breaking of con- 
tinuous translational and global gauge symmetry. The 
first mentioning of such a state goes back to Gross [2], 
who discovered a density-modulated superfluid emerging 
from a simple mean-field model for solid Helium. While 
this picture turned out too simplistic for real crystals, it 
was later shown ^ |4j that supersolidity can still exist 
under two key assumptions: (i) that the ground state of 
a bosonic crystal contains defects such as vacancies and 
interstitials, and (ii) that these defects are delocalized, 
thereby, giving rise to superfluidity. Forty years after its 
conjecture by Andreev and Lifshitz [3] and Chester [4] 
(ALC) the precise physical nature and mere existence of 
this enigmatic phase of matter has since remained under 
active debate [5]. 

In 2004 Kim and Chan provided first suggestive evi- 
dence for superfluidity in solid ^He in torsional oscilla- 
tor experiments [6 . Their discovery has sparked a host 
of new experimental activity [7 that qualitatively con- 
firmed this first observation, but also led to contrasting 
interpretations. Theoretical work has established that 
crystal incommensurability is a necessary condition for 
superfluidity [8 and that zero-point defects in ground 
state solid He are prevented by a large activation en- 
ergy [9l[T0]. In addition, Refs. piQHT? have shown that 
point-like defects experience an effective attraction that 
results in defect-clustering and phase separation, ruling 
out the possibility of defect-induced supersolidity as 
in the ALC scenario. Very recent experiments have fi- 
nally shown that the original observations [6 were caused 
by shear modulus stiffening of bulk solid He, and found 
no signature of superfluidity upon avoiding this undesired 
effect [T4 . As a result, there now seems to be consistent 
experimental and theoretical evidence for the absence of 
the long-sought supersolid phase in He. The mere exis- 
tence of free-space supersolidity in any type of quantum 
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FIG. 1: (color online) Zero temperature- limit phase dia- 
gram of two-dimensional soft-core Bosons described by eq.([T]) 
displaying the emergence of superfluid (SF) and different 
solid (NS) and supersolid (SS) phases for varying interac- 
tion strength U and density p. For the left y-d^xis the den- 
sity has been scaled by the soft-core radius Rc, while the 
right axis gives the density in units of the inverse area, 
A = v^(l -6/^0)^/2, of the unit cell of the high-density solid 
phase, corresponding to the lattice site occupation N/Ns for a 
given number of particles and lattice sites, N and iVg, respec- 
tively. The grey region labeled as NS for A > 1.5 corresponds 
to a cluster crystal. Supersolid phases with different occupa- 
tion numbers are found between two hyperbolas, defined by 
RlpU — const (dotted lines). At high densities [Ap > 3.5) 
they can be understood in terms of density modulated super- 
fluids, while superfluidity within the low-density supersolid 
lobes emerges from delocalized zero-point defects. 



solid thus remains an open question. 

In this Letter, we show that a quantum gas of bosonic 
particles interacting via soft-core potentials provides a 
prototype system that realizes all of the described phases. 
Using exact numerical techniques, we determine the com- 
plete zero-temperature phase diagram which, in addi- 
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tion to a superfiuid and normal crystal phase, contains 
a density-wave supersolidity. Most importantly, this sys- 
tem is shown to form incommensurate superfiuid ground- 
state solids — providing the first example of defect- 
induced supersolidity as conjectured by ALC [3", 'T. 

We consider a two-dimensional ensemble of N Bosons 
with density p, interacting via a pair potential of the type 
V{r) = Vo/{r^ ^R2). This interaction approaches a con- 
stant value Vq as the inter-particle distance, r, decreases 
below a critical soft-core distance i?c, and drops to zero 
for r > Rc. The limiting case 7 ^ oo yields the soft-disc 
model [15] , while 7 = 8 and 7 = 6 correspond to soft-core 
dipole-dipole [I6l [17] and van der Waals [18] interactions 
that can be realized with ultracold atoms [18l [19] or po- 
lar molecules [20l[2T]. Here, we focus on the latter case 
(7 = 6), for which the Hamiltonian reads 

where the units of length and energy are Rc and /mR^^ 
respectively, and m denotes the particle mass. In these 
units, the system behavior is controlled by the dimen- 
sionless interaction strength U = mVo/{Ji?R^) and the 
dimensionless density R^p- 

Particles with soft-core interactions have been studied 
previously in the field of soft condensed matter physics 
[22l[23], in the classical high-temperature regime. One 
of the main findings has been that pair potentials with 
a negative Fourier component [22] favor the formation of 
particle clusters, which in turn can crystallize to form a 
so-called cluster-crystal. In the quantum regime, theoret- 
ical work has so far focused on the regime of weak inter- 
actions and high particle densities [I6l [TSl [24H27] , which 
was shown to be well described by mean-field calcula- 
tions [28l[29]. In this limit, one finds strongly modulated 
superfiuid states [TS] [161 HI] 5 corresponding to a density- 
wave supersolid as originally conjectured by Gross [2 . 

Here, we investigate the so far unexplored domain 
of large interaction strengths and low particle densities, 
where strong correlations and quantum fluctuations are 
expected to become important. We employ path integral 
Monte Carlo (PIMC) simulations [30] based on the con- 
tinuous space worm algorithm [31] to determine the equi- 
librium properties of the Hamiltonian ([T]) in the canonical 
ensemble, i.e. at fixed density and temperature. An ac- 
curate determination of the ground state in the thermo- 
dynamic limit is obtained by decreasing the temperature 
to very small values sufficiently below /mR^ and by 
increasing the system size until we reach convergence of 
relevant physical observables, such as the superfiuid frac- 
tion fs and density-density correlation functions ^2(^)5 
see e.g. [T1[T6]. These quantities are used to distinguish, 
e.g., between the superfiuid phase (finite value of /g), the 
crystalline phase (large density modulations in ^2(^) and 



vanishing fs) and supersolid phase (large density modu- 
lations in g2{r) and finite /g). 

The obtained phase diagram, shown in Figjl] reveals a 
rich spectrum of phases with varying interaction strength 
and density. At small densities R^p < 0.5 we find two 
phases: a superfiuid and an insulating triangular crys- 
tal composed of = N singly occupied sites, being Ns 
and N the number of sites and the total number of par- 
ticles, respectively. The observed lobe structure of this 
crystalline region is readily understood by noticing that 
at very low densities, i.e. large inter-particle distances 
f = l/y^7fp > i?c7 the physics is dominated by the long- 
range tail of the interaction potential, V ~ 1/r^. For 
a fixed interaction strength [/ > 35, we thus find a first 
order liquid-solid quantum phase transition with increas- 
ing Ap, consistent with previous work on Bosons with 
power-law interactions [2Q [ [32H34] . However, in marked 
contrast to such systems, here the repulsion weakens 
as the particle distance approaches the soft-core radius 
Rc and declines dramatically beyond the turning point 
Ro = (5/7) ^/^i?c- Consequently, the crystal melts again 
for increasing density (f < Rq) leading to a reentrant 
superfiuid. Indeed, the resulting insulating lobe is cen- 
tered around the density at which f = R^^ as indicated 
in Fig{l] 

A distinctive consequence of the soft-core interaction 
is that the energy cost for forming close particle pairs is 
bound by Vq. This potentially enables the formation of 
crystalline phases with N > above a critical density 
where doubly occupied lattice sites become energetically 
favorable upon increasing the lattice constant. As ex- 
pected for a triangular crystal, the lattice constant de- 
creases as a = (a/3p/2)~^/^ at small densities. However, 
around a ^ lARc it increases again and settles to a 
density-independent value of ao — 1.6 Rc upon further 
increase of p. The corresponding volume of the unit cell 
A = Y^3/4ao provides a measure of the lattice occupancy 
N/Ns = Ap, which is shown in Figjl] The transition to 
cluster crystals occurs at Ap ^ 1.5, which indeed coin- 
cides with the critical density for crystallization of the 
reentrant superfiuid phase. 

Around this density, a thin region of phase separation 
is found to lay in between the cluster crystal and the su- 
perfiuid phase. Figj2|b) shows a typical example for the 
particles density in this region. Two distinct coexisting 
phases can be recognized: a crystal phase with exactly 
two particles per site (upper part of the figure), and a 
superfiuid phase (lower part). We have carefully checked 
that the occurrence of this phase-separated state is not 
an artifact of the simulations, by performing accurate an- 
nealing and by choosing different initial conditions, such 
as random and different crystalline configurations. 

Above incommensurate lattice occupations N/N^ ^ 
1.5, the direct liquid-solid quantum phase transition is 
replaced by a first order transition from a superfiuid 
to a supersolid phase. The SS phase is approximately 
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FIG. 2: (color online) (a)-(f) Two-point correlation function g2{r) (angle- averaged) and density distribution p(r) for pR^U — 32 
and different values of Ap^ indicated in panel (g), showing the superfluid fraction fs as a function of Ap across the different 
phases of Figjl] The horizontal dotted line shows the p ^ oo limit of fs for a density- wave supersolid, obtained from a mean 
field evaluation according to Ref. ^35j . 



found to occur between the two hyperbola defined by 
R^pU = a, with a ~ 28 and a ~ 38, respectively (see 
dotted lines in Figjl]). These two lines are derived from 
the weak interaction limit {U and p ^ oo with 
a = const.), where mean-field theory predicts a transi- 
tion to a density- wave supersolid that is determined only 
by the value of a [18 . While mean-field predictions are 
essentially exact for densities Ap > 4 (see Figjl]), the sit- 
uation is dramatically different at smaller densities where 
the discrete nature of the particles plays a significant role. 
This gives rise to the emergence of supersolid regions with 
a lobe structure, which vanish at commensurate lattice 
occupations N/N^ = 2 and N/N^ = 3. There, we find a 
direct transition between a superfluid and an insulating 
solid phase. 

This behavior is illustrated in Figj2j where the super- 
fluid fraction, /g, is shown as a function of Ap for a = 32, 
i.e. in between the two dotted lines in Figjl] As shown 
in Figj2]^c) and (e), at Ap — 2 and Ap — 3 one finds a 
commensurate crystal with exactly N/N^ = 2 [Figj2]^c)] 
and N/N, = 3 [Fig§e)] per lattice site, respectively, and 
vanishing superfluidity [Figj2]^g)]. Importantly, the crys- 
tal structure in between these two densities is practically 
unchanged, as seen by comparing the particle density and 



density-density correlation function g2{r) (averaged over 
angles) in Figsj2]^c)-(e). However, the incommensurate 
lattice fllling N/N^ and the resulting fluctuations of in- 
dividual occupations enables particles to tunnel between 
the sites. This gives rise to a non- vanishing superfluid 
fraction of the crystal, which can assume sizable values 
of /s = 0.3 for N/N^ ^ Ap = 2.5. 

For higher densities Ap > 3, the scenario described 
above changes considerably. As shown in Figj2]^g) 
the superfluid fraction approaches a constant, density- 
independent value /s ~ 0.24 with increasing p. In par- 
ticular, for an average commensurate fllling N/N^ = 4 
there is no direct phase transition between a superfluid 
and a solid insulating phase, and instead the supersolid 
phase persists with no signiflcant difference to the case of 
incommensurate lattice occupancies N/N^ ^ 4. This be- 
havior signals a crossover to the regime where the super- 
solid phase can be understood in terms of a density mod- 
ulated superfluid [2j [151 [161 [181 [28] , where the discrete na- 
ture of the particles becomes irrelevant. In this limit, the 
superfluid-supersolid quantum phase transition is well 
captured within a mean-fleld description [151 EHl [29] , 
which predicts a transition point at a = R^pU = 28.2 
as well as a superfluid fraction /g = 0.23 that is solely 
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defect fraction, /def 

FIG. 3: (color online) Superfluid fraction in a doubly occu- 
pied solid (N/Ns = 2) as a function of the fraction of defects 
in the form of singly occupied sites. Here we fixed [/ = 31, 
A^s = 72, and = 139 ^ 144. The inset shows the particle 
density (color code) for fdei = 0.03, obtained from a PIMC 
configuration along with the particle positions for a single 
imaginary-time slice. There are two singly occupied sites, 
which are highly mobile as indicated by the corresponding 
imaginary-time trajectories. 



determined by a. As shown in FigsjT] and [2|g) this is 
well confirmed by our PIMC results for Ap ^3.5, sug- 
gesting that the transition to density-wave supersolidity 
takes place at a small particle number per lattice site. 

The most interesting behavior takes place around the 
superfluid-solid quantum phase transition at N/N^ = 2. 
Figure |3] provides a more detailed look at the transition 
between the insulating crystal and the supersolid phase 
close to R^pU ^ 28 {U = 31). Starting from the insu- 
lating solid with doubly occupied lattice sites, we succes- 
sively remove a small number of particles from randomly 
chosen sites and monitor the superfluid fraction of the 
resulting new groundstate obtained from our PIMC sim- 
ulations. Removing a small number of particles does not 
cause structural changes of the ground state but rather 
creates a small fraction /def = {2Ns — N)/Ns of zero-point 
crystal defects in the form of singly occupied sites. An 
analysis of our PIMC configurations shows that defects 
do not cluster and instead delocalise, as illustrated in the 
inset of Fig|3] Indeed, we find a finite superfluid fraction 
even for small defect concentrations, which increases lin- 
early with /def. We have checked that this behaviour 
is pertinent to the groundstate and not to a metastable 
configuration by running simulations with different initial 
conditions, including clustered defects. This behaviour 
is thus consistent with defect-induced supersolidity ac- 
cording to the ALC scenario, and constitutes the central 
result of this work. 

Supersolidity in this system is the consequence of two 
unique features of soft bosons, (i) The energy cost for 



forming close particle pairs is bound by Vq, which fa- 
cilitates the formation of cluster crystals that naturally 
entail zero-point defects, (ii) The dynamics and interac- 
tion of these defects differs fundamentally from those of 
conventional solids. In the latter case, vacancies and in- 
terstitials induce displacement fields that lead to purely 
attractive defect-interactions [10i|T2j|13|[36] and, thereby, 
prevent a delocalization of defects [10]. In cluster solids, 
on the other hand, defect interactions are purely repul- 
sive, since two interstitials interact via essentially the 
same underlying particle interaction V(r) of eq. ([T]). The 
same applies to vacancies due to particle-hole symmetry 
around integer N/N^. In the present case, the transition 
between these two regimes is controlled by the particle 
density. For Ap > 1.5 delocalized zero-point defects al- 
low for the formation of supersolid phases. Below this 
density particles do not explore the soft core part of the 
interaction potential, such that defects are attractive and 
supersolidity is absent consistent with the results of Ref. 
[10 . Around the transition region Ap ^1.5 neither pic- 
ture applies, and one observes separation between a su- 
perfluid and a doubly occupied, insulating cluster solid. 
Preliminary calculations based on path integral Langevin 
dynamics [37] suggest that in this region structural and 
dynamical heterogeneity gives rise to a quantum glass 
phase at finite temperature. 

Having identified a physical system that facilitates 
defect-induced supersolidity, we hope that this work 
will provide useful guidance for future experiments and 
initiate further theoretical explorations. An important 
question concerns the general features of the interaction 
potential that are required to maintain the type of 
supersolid states described in this work. While the 
emergence of density-wave supersolids is largely insen- 
sitive to the detailed shape of the soft-core interaction, 
the low-density physics described in this work may 
be strongly affected. In fact, it seems reasonable to 
expect an interesting competition between intra- and 
inter-site interactions within the self- assembled crystal 
that will strongly depend on the long-range tail of the 
particle interactions. While we have focussed here on 
two-dimensional systems, the role of the dimensionality 
represents another outstanding issue, and in particular 
its role for frustration effects that may be decisive for 
defect delocalization. While the considered interactions 
do not straightforwardly occur in natural crystals, we 
note that so-called Rydberg dressing [H [TQI [38H4Q] 
of atomic Bose-Einstein condensates constitutes a 
promising approach for an experimental realization. 
Ultimately, the characterization of a broad class of 
suitable interactions and conditions will greatly ease 
the search for viable implementations. At this point, 
ultracold atoms and molecules appear most promising, 
as they offer a great deal of flexibility in designing 
various types of long-range interactions |20l . 
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